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ABSTRACT 


Zinc is a major residue impurity in the preparation of solar grade 
silicon ma* rial by the zinc vapor reduction of silicon tetrachloride. 
This paper projects that in order to get a 17% AMI cell efficiency for 
the Block IV module of the Low-Cost Solar Array Project*, the concen- 
tration of the zinc recombination centers in the base region of silicon 

11 3 

solar cells must be less than 4x10 Zn/cm in the p-base n+/p/p+ cell 
11 3 

and 7x10 Zn/cm in the n -base p+/n/n+ cell for a base dopant impurity 

14 3 

concentration of 5x10 atoms/cm . If the base dopant impurity concen- 

15 3 

tration Is increased by a factor of 10 to 5x10 atoms/cm , then the 
maximum allowable zinc concentration is increased by a factor of about 
two for a 17% AMI efficiency. The thermal equilibrium electron and hole 
recombination and generation rates at the double acceptor zinc centers 
are obtained from previous high field measurement;; as well as new 
measurements at zero field described in this pc.; r . These rates are 
used in the exact d.c. circuit model to compute the projections. 


*Seo for example JPL/DOF. 8101-14 1 p.21, a report distributed at the 
14th Photovoltaic Specialist Conference, January 19S0. 



I. INTRODUCTION 


The objective of this program i3 to determine the effects of 
impurities and defects on the performance and permanence of silicon solar 
cells. It includes theoretical (computer model) and experimental studies 
of the effects of impurities on the properties of silicon intentionally 
doped with specific impurity elements, as well as the effects of these 
impurities on the impurity related energy level positions, the concentration 
of these energy levels and the recombination-generation properties of electrons 
and holes at these energy levels. 

This third technical report contains a theoretical projection of 
the maximum level of zinc recombination centers in silicon which is allowable 
to give a AMI silicon solar cell efficiency of 17%. 

The thermal recombination and generation rates of electrons and holes 
at the double acceptor zinc centers in silicon are obtained from published 
results in the literature as well as new measurements described in this report. 
These values, considered to be the currently best values at 300K, are used 
in the transmission line equivalent circuit model to compute the performance 
of silicon p+/n/n+ and n+/p/p+ cells doped with zinc from 10^ to 10^ Zn/cm^. 

Chapter II gives a detailed description of the zinc recombination 
and generation rate data, the material and device parameters used in the 
theoretical computation, and the performance of the silicon solar cells 
computed using the one-dimensional exact equivalent circuit model and these 
parameters. The AMI solar spectral density, the silicon absorption coefficient 
and silicon surface reflection coefficient are tabulated in the appendix. 

The silicon properties, the intrinsic carrier concentration and the electron 
and hole mobilities, used in the theoretical computation are also given in the 


appendix. 



II. EFFECT OF ZINC IMPURITY ON SILICON SOLAR CELL EFFICIENCY 


1. INTRODUCTION 

In producing low cost silicon stock by zinc vapor reduction of silicon 
tetrachloride in a fluidized bed, a large amount of zinc impurity remains in 
the granular silicon product [1], The residue zinc will give high concentrations 
of zinc recombination centers in single crystals or sheets grown using this zinc 
reduced silicon stock. This analysis provides a guideline on the maximum 
amount of zinc allowable in a silicon solar cell fabricated frcm such a single 
silicon crystal or sheet for a desired AMI energy conversion efficiency. 

2. RECOMBINATION PROPERTIES OF ZINC CENTERS IN SILICON 

Zinc is a double acceptor in silicon [2,3]. The first acceptor level is 
located at 326 meV above the silicon valence band edge and the second acceptor 
level is located at 664 meV above the silicon valence band edge. These were 
determined by Herman [4,5] using the transient capacitance and current methods 
[6,7,8], Strong and nearly equal-amplitude voltage stimulated capacitance 
transient (VSCAP) [8] and thermally-stimulated capacitance (TSCAP) [9] signals 
are observed due to thermal emission of holes trapped at these two levels. The 
VSCAP peaks are illustrated in Figure 1 

Detailed and extensive measurements of the thermal capture and emission 
rates of electrons and holes at these zinc levels have been made by Herman [4,5] 
in diffused silicon p+/n and n+/p diodes which were also zinc diffused in a 
closed quartz ampoule, however, these measurements were made at high electric 
fields in reverse biased p/n junctions. In order to predict the performance of 
silicon solar cells containing zinc, recombination rates at low or zero electric 
fields are needed. Some of the high field data can he extrapolated to low 
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Capacitance Transient Signal, C 9 -C 1 , (100 mpF) 



Figure 1 Computer plot of sampled voltage stimulated capacitance (VSCAP) 

as a function of diode temperature of a zinc-diffused silicon 

n+/p diffused diode. The two peaks are from thermal emission 

of holes trapped at the two zinc acceptor levels, E +326 and 

F. +664 meV above the valence band edge of silicon, 
v 
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fields but Additional low field data are still needed to give sufficiently 
accurate thermal capture rates of electrons and holes at the zinc levels at 
low fields in order to predict the solar cell performance. These were made 
and the results of the thermal capture and emission rates at the two zinc 
acceptor levels are listed in Table 1 and compared with those values used 
in a recent study of the dark forward current-voltage (I-V) characteristics 
of zinc doped silicon diodes. [10] The methods by which these new values were 
obtained are discussed In the following two subsections. 

2.1 FIRST ACCEPTOR LEVEL (E y + 326 meV) 

The thermal recombination kinetics at this level are characterized by 
the two capture rates, Cp^ and c n Q, and the two emission rates, and e^. 

The numeral subscript, 0 or 1, Indicates the number of electrons trapped at the 
zinc center pfu'oA. to the Indicated electron or hole transition [11]. The 
reference state, subscript 0, Is the electrically neutral charge state. Thus, 
c^ is the capture rate of holes by a negatively charged zinc center. The 
determination of their values at room temperature, 23. 4C or 296. 6K, is described 
for each of the rate constant given below. The temperature 296. 6K is used here 
since It corresponds to a convenient value of the intrinsic carrier density of 
n^*l. 000x10^ carriers/cm^ . 



The thermal hole capture rates at singly negatively charged zinc center, 

c ,, were measured by Herman [5] at 83K as a function of the electric field from 
Pi 

4 5 

2.5x10 to 10 V/cm. This is extrapolated to give a zero field value of 
-7 3 

(4 t 2)xl0 cm /sec. Zero field capture rates tree also measured by 
us using the voltage stimulated capacitance transient (VSCAP) [8] at 112.57 
and 122. 7K with values of 1.6x10 ^ and 3.9x10 ^ em^/sec respectively. 


S- 


Table 1 


The Thermal Ealaalon and Capture Rate* at the Double 
Acceptor Zinc Centers in Silicon at Room Temperature 


Energy Level ■ Ey + 664 tneV 


Rates 

This Work 
296. 6K 

Dark T-V* 
297. OK 

Sources 

3 

c nl (cm ft) 

-12 

2 x 10 

5xl0~ U 

Exp. 

*n 2 a/ *> 

0.41 

11 

Computed 

3 

c » (cm ft) 
P* 

6 xl 0' 8 

6 xl 0 -8 

Exp. 

* p i (1/ *> 

29 

29 

Exp. 


Energy Level « Ey + 326 meV 


Rates 

This Work 

Dark I-V® 

Sources 


296. 6K 

29 7. OK 


c „ 0 (cm 3 /s) 

2 xl 0' 8 

-9 

2 x 10 

Exp. 

*«1 (1/s) 

5xlO ~ 3 

9xlO ' 3 

Computed 

c (cm 3 /s) 

P 1 

5xl0 " 8 

9xl0 -8 

Exp. 

e p 0 (1/,) 

2 xl 0 7 

2 xl 0 7 

Exp. 


Rate data from reference 10. 


- 6 - 








































In our recent study of the dark forward I-V characteristics of zinc 

diffused silicon diodes [10], it was assumed that c . was independent of temperature 

\ P 

so that the 83K value of Herman just quoted was used at room temperatures. 

Although the new values obtained from direct capture rate measurements at 123K 

showed that the capture rates are nearly constant from 83 to 123K, we prefer 

an extrapolation to 296. 6K be made by following the temperature dependence of 

the c^ capture rate since in principle we expect C p^ <C p 2 * ^ C pl were taken 

to be constant, then at 296K, we would have c ^ > c^ while the experimental 

data at T< 125K showed c , <c The sources of c . data are described in the 

pi p2 p2 

next subsection. 

-8 3 

This extrapolation gives a value of Cp^=5xl0 cm /sec at 296. 6K which is 

-8 3 

two (2) times smaller than the previous value of 9x10 cm /sec [10]. 

(B) Thermal Electron Capture Rate , c^q 

The thermal electron capture rates at the neutrai zinc center, c n Q» were 
measured by us recently using the diode forward-injection-reverse-recovery 
method [12] since the VSCAP experiment cannot set the initial charge states of 
the zinc centers to the neutral charge state to observe this electron transition. 

The use of this diode current recovery method requires a careful interpretation 
of the result based on the recombination kinetics of a double acceptor level. 

The necessary theoretical model was given previously [13]. The high level 

2 

lifetime at the experimental forward injection current of 10 mA or 1200 mA/cm 
had a nearly constant value of 633 ns in a zinc diffused p+/n diode and 643 ns 
in a zinc diffused n+/p diode. The diodes used in these experiments were those 
fabricated by Herman [4]. The nearly constant and equal high level lifetimes 
in the n-type and p-type bases are consistent with the theoretical considerations 
which lead to the conclusion that electron capture by the neutral zinc center 
should be the dominant recombination transition which controls the high level 


- 7 - 



lifetine in the base region. This is arrived at as follows. From the two-level 
theory [13], the high level steady-state lifetine or transient lifetime is given 
by 

t H * <c nO c nl + c nO c p2 + VV' [C nO c p2 (c pl + 'nl^TT 1 (1) 

In general, we would expect c n i <<c pi and c n i <<c p2 s * nce c n i is associated with 
a repulsive potential from the singly negatively charged zinc center while 
both c^ and c^ are associated with attractive potentials from singly and 
doubly negatively charged zinc centers. The experimental data of c ^ described 
in the next subsection indeed confirms these inequalities. Thus, the high level 
lifetime given by Equation (1) simplifies to 

t H ~ (c n0 + c pl )/c nO c pl = T n0 + T pl <2) 

Mhere T n<T 1/c nO N TT and T pl =1/c pl^TT’ 

The simplified result given by Equation (2) shows that the high level 
recombination lifetime in a quasi-neutral base region is nearly entirely 
determined by the recombination through the lower or first acceptor level. 

This result is quite general and it also applies to a double donor center such 
as sulfur in silicon. However, it is generally not applicable to a two-level 
amphoteric impurity such as gold which has a donor and an acceptor level in 
silicon, except for special circumstances. For the amphoteric center, we have 
no assurance that one of the capture rates is substantially smaller than the 
other three, however, the case of gold is close to this situation, making the 
lower gold donor level the dominant recombination level at high injection levels. 

For the first zinc acceptor level considered in this section, the 

1 A 3 

zinc concentration in the diodes is N TT »10 Zn/cm whose maximum value may be as 
high as 1.7xlO x Zn/cni at the junction due to pile up and enhanced solubility 
in the phosphorus diffused n+ emitter layer. Using the smaller and bulk value, 
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and c pl “5xlO cm /sec just obtained in section (A), the electron capture rate 
at the neutral zinc center is 


C n0 ’ 1/T nO N TT * 1/(t H “ T pl )N TT 


1/(643-200)x10" 9 x10 14 - 2.25xl0“ 8 cm 3 /sec 


(3) 


-9 3 

This is ten ( >) times larger than the previously used value of 2x10 cm /sec 
[10]. It would increase the efficiency reduction in n-based silicon solar cell 
by the presence of zinc. 

(C) Thermal Hole Emission Rate , e p Q 

The thermal hole mission rate, e^, is computed from the experimental 

Arrhenius equation obtained by Herman at low temperatures [4], 

e n - 7. 4xl0 12 (T/ 300) _1 exp(-326/kT) (4) 

pU 

where k«0. 08617087 meV/K. Herman found a large electric field dependence of 

this hole thermal emission rate. The T * dependence of the pre-exponential or 

-2.5 -2.5 1.5 

collision frequency factor is selected since c^ « T so c^Ny “ T T 

-1 -2 
* T . A different temperature exponent, such as T , will change e^ by less 

than a factor of 2. However, the magnitude of e n , if varied by a factor of 2 

pO 

or even more, will have little effect on the predicted solar cell performance 
since it depends mainly on the capture rates and not the emission rates. The 
value computed from Equation (4) is 2.2xl0 2 1/sec at 297K and about 2xl0 2 1/sec 
at 296. 6K which is used in this work. 

(D) Thermal Electron Emission Rate , e^ 

The thermal electron emission rate from the negatively charged zinc center 

is computed from the mass action law at thermal equilibrium since under forward 


'nl 


bias, the quasi-neutral base region of a solar cell is nearly at thermal 


10 -3 


equilibrium. The mass action law is e n i e pQ^ C n0 C pl " n i‘ ^ s * n 8 n ^*1.0x10 cm 


-3 

at 296. 6K, we have e^ = 5x10 1/sec. 
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2.2 SECOND ACCEPTOR LEVEL (Ey + 664 maV) 

The thermal recombination kinetics at this level are characterized by 
the two capture rates, c^ *nd c^, and the two emission rates, and e^. 

(A) Thermal Hole Capture Rate . c p2 

“8 3 

A low field value of 6x10 cm /sec is extrapolated from Herman's high 
field data at three temperatures [5], 1.4xl0 -6 at 135K, 7. 5x10“ 7 at 159K, and 
3.0x10 7 at 190K. This is the same extrapolated value used by us previously 
in diode dark current-voltage characteristics studies [10]. 

(B) Thermal Electron Capture Rate , c^ 

This electron capture rate is expected to be very small since the electron 
is captured into a repulsive Impurity potential of the negatively charged zinc 
center. The previous values were obtained by Herman [5] at low temperatures and 
in high electric fields, (1-8) xlO* V/cm. His values were 4.6x10 33 cm 3 / sec at 
138K, 4.4x10 33 cra 3 /sec at 150K, 5.2x10 33 cm 3 /sec at 165K and l.lxlO -33 cm 3 /sec 
at 181K. Since the capture is by a repulsive potential well, we would expect 
the capture rate to be highly dependent on the electric field. At lower fields, 
the impurity potential barrier would not be lowered as much and the capture 
rate should be substantially reduced from the high field value. Thus, new 
measurements of the thermal equilibrium capture rates were made by us 
using the VSCAP method [8]. Data were obtained at three temperatures, giving 
4x10 33 cm 3 /sec at 230K, 6.2x10 33 cm 3 /sec at 250K and 1.1x10 3 ^ cm 3 /sec 
at 270K. These data can be fitted to the Arrhenius equation since the repulsive 
potential well gives rise to a thermally activated process. The result is 

C nl " (3 * 7 ± 2,8 ) xl0 ” 10 exp[-(136±16)/kT] cm 3 / sec (5) 

-12 3 

Extrapolating to 296. 6K gives a value of 2x10 cm /sec. 

This thermal equilibrium value is twenty-five (25) times smaller than the . 

-11 3 

high field value at low temperatures previously used, 5x10 cm /sec. The large 
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high field value at low temperatures is probably due to the lowering of the 
repulsive potential barrier by the high electric field. At low fields, the 
barrier height is about 136 meV as deduced by the thermal activation energy 
of the electron capture rate given by Equation (5) . 

(C) Thermal Hole Emission Rate, e 

Pi 

The thermal hole emission rate from the singly negatively charged sine 


center, e^, is again obtained from the experimental Arrhenius equation 

obtained by Herman at low temperatures [4], 

e pl « 5.36 x 10 12 (T/300)” 2 exp(-664/kT) (6) 

4 5 

who found essentially no electric field dependences from 10 to 10 V/cm and 

-2 

in the range of 200 to 360K. The T depenaence is selected here since 


4 

:T , however, a different power will change e p ^ 




which would have no effect on the recombination rate under solar cell operating 


conditions. At 297K, Equation (6) gives e p j”28.94 1/sec. If the T~* dependence 

is used with AE-641 meV, then e ,-28.64 1/sec which differs very little from 

Pi 

28.94 given by Equation (6). 


(D) Thermal Electron Emission Rate , e^ 

The thermal electron emission rate at the doubly negatively charged 

2 

zinc center is computed from the mass action law, e n 2 e p i “ c n i c p 2 n i 
gives e^ “ 0.41 sec 
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THEORETICAL AND COMPUTER MODEL 


Hie performance of silicon solar cells was computed by numerically solving 
the six Shockley Equations for two-level recombination centers [14]. A 198 lump 
section transmission line equivalent circuit model was synthesized from these 
equations [13,15]. A photocurrent generator was added to each of the lump 
sections to represent the solar Illumination. Numerical solutions of the quasi- 
Fermi potentials for electrons and holes and the electrostatic potential at each 
of the 198 sections were obtained by standard circuit analysis technique involving 
the diagonlzation of a 198x3 sparse matrix. Infinite interface recombination at 
the front and back surfaces were simulated by short-circuiting the three potential 
nodes at each of the two interfaces. The solution procedure began by assuming 
a distribution of the electrostatic potential at equilibrium (zero applied voltage 
in the dark) where the quasi-Fermi potentials have a constant value. The errors 
of thi3 initial guess of the electrostatic potential were computed from the small- 
error transmission line circuit model and added to the initial guess to get a new 
set of 198 electrostatic potentials. This iteration procedure was continued 
until a specified accuracy was reached by all 198 electrostatic potentials. 

In the next step, a small photocurrent generator was added between the 
electron and hole quasi-Fermi potential nodes at each of the 198 positions while 
the applied voltage was kept at zero. The converged dark solution of the 
electrostatic potentials just obtained were used as the initial guess. The 
interation procedure was then repeated to get the converged solutions of the three 
potentials at each of the 198 positions. The photocurrent generator was then 
increased and the iteration calculation repeated until the illumination reached 
the AMI intensity, giving the AMI short-circuit solutions. Then, the applied 
voltage was stepped in 50 mV increments and the iteration repeated until an applied 
voltage where the computed diode current changed sign. At this point, the 
solution had passed the or n-circuit condition and the iteration was stopped. 
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The accurate ma ximu m power efficiency and open circuit voltage vara 
o bt ai ne d by a fifth order polynomial extrapolation of the numerical eolutlone 
Juat obtained in the vicinity of theae two voltages. The extrapolations were 
made twice* each tine with a new set of quasi-Fermi potentials and electrostatic 

potentials computed at the extrapolated maximum power voltage or open circuit 
voltage. 

An illustration of the numerical procedure was given by Maes and Sah [16] 
who computed the d.c. dark I-V characteristics of gold diffused silicon diodes 
and by Chan and Sah who computed the d.c. dark I-V characteristics of sine 
diffused silicon diodes recently [10]. The differences between these dark 
forward I-V calculations and the illuminated forward I-V calculations are 
the addition of the photocurrent generator and the I-V locus used in the 
calculation. The illuminated I-V locus of calculation was just described. 

To compute the efficiency of silicon solar cells containing zinc* the 
following diode material and device properties were selected. 

(1) Zinc concentration is assumed constant over the diode. This 
assumption may differ from practical situations where the high concentration 
of the phosphorus or boron diffusant in the diffused emitter or the back 
surface field layer may enhance or depress the solubility of zinc. The 
enhanced solubility will only slightly decrease the cell performance since 
the results, to be presented, show that base recombination dominates the cell 
performance for practical base lifetimes for 10^ to 10^ Zn/cm^. 

(2) Interband Auger recombination was included in the computation to 
give a more realistic situation when the zinc concentration was extremely low. 
However, no other base recombination other than zinc is Included. 

(3) Both the n+/p/p+ and p+/n/n+ back surface field cell structures 
are analyzed numerically. The thickness of the cell is taken as 250 pro. The 
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diffused emitter junction le located at 0.25 urn from the front aurface and 

the diffused back aurface field layer concentration dropa to bulk dopant value 

at 1.0 pm from the back surface. Measured boron and diffusion profiles are 

fitted analytically and used here [17,18]. The phosphorus profile is fitted to 

exp(-X^) due to its rather flat plateau and sharp drop off. The boron profile 

2/3 

is fitted to (1 - X ) again simulating a flat plateau and sharp drop off. 

The rather abrupt phosphorus profile increases emitter recombination and 

2 

reduces the short circuit current by nearly 2 mA/cm while the boron emitter 

2 

profile reduces the short circuit current only about 0.2 mA/cm . Neither of 
these profiles have significant effects on the open circuit voltage which is 
determined mainly by recombination at the zinc centers in the base region. The 
profiles employed are given below. 

Naa - Ndd - C 0 fl - x’%] - C B - C L exp[-(L-x) 6 /L 2 ] (p+/n/n+) 

Njjjj - N aa - ^expl-x 6 /^] - C B - C L U - (L-x) 2/3 /L 2 ] (n+/p/p+) 

where L-250 pm and and L 2 are determined by the Junction depth of 0.25 pm 

20 -3 

for the emitter and 1.0 pm for the back surface field. A value of 2.5x10 cm 

is used for both the front and back surface concentrations, Cq and C^. The 

bulk concentration is taken as 5x10** cm” 3 . For p-base or n+/p/p+ cells, 

this corresponds to a base resistivity of 25 ohm-cm. For n-base or p+/n/n+ 

cells, this corresponds to a base resistivity of 10 ohm-cm. An increase of 

15 -3 

the base doping concentration by 10 to 5x10 cm will increase the AMI 
maximum efficiency by about 1% (from 16% to 17%). 

(4) Antireflection coating is assumed to eliminate surface reflection 
completely. For a bare cell surface without antireflection coating, the 
snort circuit current is reduced by about 36%. 

(5) The AMI solar spectral density tabulated by Thekaehara [19] at the 
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ambient condition of 20 an H^O and 3.4 on 0^ with a total Incident power of 
2 

88.92 mW/cm ie employed. This la aomewhat amaller than the industrial 

2 

standard of 100 mW/cm . The photocurrent generator located at position 
for the lump section located between x^ and x^ Is given by 
J 0 (x l ) " q j ^ <^ hc > p <°* X >I«" C “ l ' •“® t *]dX/(x 2 -x 1 ) 

2 

where P(0,X) is the incident AMI solar power spectral density (W/m y), X the 
photon wavelength* h the Plank constant, c the velocity of light in air, X n 
the cut-off wavelength of silicon determined by the silicon absorption 
coefficient, a, and approximately given by X^-hc/F^, *nd Xo the cut-off 
wave-length of the solar spectrum. This photocurrent generator gives the 
value averaged over the lump section whose thickness is (x^-x^). Such an 
average insures that computed short circuit current does not exceed the 
available photocurrent due to numerical inaccuracy. 

(6) The silicon absorption coefficient as well as the reflection 
coefficient described in (4) were obtained by interpolation of the data 

given in the literature. The absorption coefficients were numerical interpolation 
of values read off enlarged figures of the data of McFarlan, et. al. [20] and 
Philipp and Taft [21]. The reflection coefficients were obtained similarly from 
the data of Philipp and Taft [21]. These are given in the table in the appendix. 
The four significant figures listed do not imply absolute accuracy but serve as 
a constant base for numerical comparison of the performance of different cell 
geometries and material properties. 

(7) The properties of the silicon material used in this analysis are also 
given in the appendix. These include the intrinsic carrier density, n^, which 
contains energy gap narrowing due to Debye screening by the high density of 
electrons or holes in the diffused layers, and the lattice (acoustical and 
optical or intervalley phonon) and impurity scattering limited mobilities. 
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PROJECTED AMI PERFORMANCE 


Th« open circuit voltage, short circuit current, fill factor and 
maximum efficiency at AMI illumination were numerically obtained. The 
fill factor la about 3.8 and varies by 0.02 among the cases and hence are not 
displayed in figures while the other factors are all given graphically. 

The maximum AMI efficiencies are shown in Figure 2. Two of the 
four curves, la and 2a, contained in ter band Auger recombination while the 
other two did not. Curves 1 and la are for the xrt7p/p+ cells and 2 and 2a 
are for the p+/n/n+ cells. 

It is evident that the p+/n/n+ cell has a slightly higher efficiency 

than the n+/p/p+ cell when the Auger recombination is Included. The small 

difference arises mainly from the interband Auger recombination in the p+ and 

n+ emitters. These results show that in order to have a cell efficiency of 

11 3 

171, the zinc concentration must be less than about 4x10 Zn/cm in the 

11 3 

p-base or n+/p/p+ cell and less than about 7x10 Zn/cm in the n-base or 
p+/n/n+ cell. 

If the base dopant impurity concentration is increased by 10 to 5x10^ 

the maximum allowable zinc concentration is increased by a factor of about 2. 

15 3 

Thus, for the p-base cell of 5x10 Boron/cm (3 ohm-cm), the maximum zinc 

11 3 

concentration for 17% efficiency is about 8x10 Zn/cm . For the n-base 
15 3 

cell of 5x10 Phosphorus/cm (1 ohm-cm), the maximum zinc concentration for 
17% efficiency is about 1.5x10*^ Zn/cm^. 

The effects of zinc concentration on the open circuit voltage and short 
circuit current are shown in Figures 3 and 4, The open circuit voltage is not 
affected by the Interband Auger recombination in the diffused emitter and 
has nearly the same values in p-base and n-base cells as shown in Figure 3. 

The Auger recombination in the emitter is not Important unless the zinc 
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Figure 3 Theoretical AMI open circuit voltage ee a function of Zn 

concentration at 296. 6K for back aurfece field silicon eolar 
cells with perfect antireflection coating. Solid curve* are 
for n+/p/p+ cells and dashed curves p+/n/n+ cells. Curves 
labeled a Include Interband Auger recombination. 
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Theoretical AMI short circuit current as a function of Zn 
concentration at 296. 6K for back surface field ail 'con solar 
calls with perfect ant iref lection coating. Solid curves are 
for n+/p/p+ calls and dashed curves p+/n/n+ calls* Curves 
labeled a Include interband Auger recombination. 




11 3 

concentration la lass than 10 Zn/cm in the boat. Whan it la greater than 

10** In the base, base recombination donlnatea the open circuit voltage. 

The abort circuit current at low sine concentration or large base 

lifetime Is mainly determined by emitter r ec omb in ation as indicated by the 

four curves In Figure 4 for the short circuit current. Znterband Auper 

recombination in the n+/p/p+ cells has a significant effect. Zt reduces the short- 

circuit current by about 2 mA/ca . This cones about from the high phosphorus 

concentration la the diffused emitter which Is nearly constant until It reaches 

the junction where It drops off abruptly. Znterband Auger recombination has a 

much mailer affect on the short circuit current In the p+/n/n+ cell, reducing 

2 

It only by about 0.2 aA/cm . Reducing the surface concentration of the diffused 

emitter will reduce the Interbend Auger recombination loss. The results shown 

In Figure 4 indicate that the short circuit current begins to decrease due to 

12 3 

base recombination when the sine concentration is higher than about 10 Zn/cm . 

5. CONCLUSIONS 

Nvaaerlcal calculations using the exact one-dimensional circuit model 

show that sine recombination centers reduce the performance of silicon solar 

cells significantly due to recombination of electrons and holes In the bass 

region of the cells. The reduction Is nearly the same for p+/n/n+ and n+/p/p+ 

cells. To get an AMI efficiency of 17X, the sine concentration In the bass 

of the n-base cell must be less than 7x10** Zn/cm® and In the p-base cell, lees 

than 4x10** Zn/cm® for a base dopant concentration of 5x10** cm - ®. Increase in 

*5 -3 

the base dopant concentration by a factor of 10 to 5x10* cm would Increase the 
maximum allowable sine concentration by a factor of two to maintain a 17% AMI 
efficiency. 
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APPENDIX 


In Table 2, the soler power spectral density, P(W/n 2 -y) , Is given as a 
function of photon wavelength, X(y), and energy E(eV). The AMI solar power 

| 

spectral irradiance is obtained by Thekaekara [19] at the atmospheric condition 
of 20 mm H^O and 3.4 mm 0^. Some numerical interpolation of Thekaekara* 8 data 
is made in preparing the table. 

The silicon absorption coefficient as a function of photon energy or 

wavelength is also tabulated in Table 2. The data are obtained by interpolation 

of the data given by MacFarlane, et. al. [20] near the threshold energy and by 

Philipp and Taft [21] in the higher energy range. Although values were read 

off enlarged pictures of the figures given by MacFarlane, et. al. and Philipp 

and Taft, the four significant figures given do not imply absolute accuracy 

but are used to provide a constant base for comparison of performance of 

different cell geometries and material. The absorption tail does provide 

2 

several tenths of mA/cm of short circuit current due to its deep penetration. 

Thus, a constant set of values is essential to provide a constant base for 
comparison. 

The surface reflection coefficient, R(X), and the transmission 
coefficient, 1 - R(X), as a function photon wavelength or energy is also 
tabulated in Table 2. These values are obtained by interpolation of the 
published data of Philipp and Taft [21] whose figures were enlarged and the 
values visually read off and then least square fitted to a polynomial. The 
four significant figures given in the table again provide a constant base for 
comparison purposes and do not imply absolute accuracy. 

In the current calculations presented in section 4, a perfect anti- 
reflection coating is assumed so that the reflection coefficient' is assumed 
to be zero or the transmission coefficient is assumed to be unity. 
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Table 2 AMI solar spectral irradiance, P, silicon absorption 

coefficient, a, and silicon surface reflection coefficient. 
R(A), as a function of photon wave length, X. 88.92 vSH/ctr* 
(part 1, 2900-6000A) 


X(U) E(«V) a(ce l ) P(W/n 2 u> 1-Rm 
. 29000+00 .42750401 .21190407 .00000400 .33290-00 
.29500400 .42030401 .10000407 .00000400 38810-00 

. 30000400 . 41330401 . US''^407 .41000401 . 38810-00 
.30500400 . 40450401 .15* T -07 .11400402 . 40650-00 
.31000400 . 40000401 .14*. . 30500402 . 41900-00 

.31500400 .39360401 . 13* 40407 .79400402 .42850-00 
.32000400 .30750401 .12700407 .20260403 .43660-00 
.32500400 . 30130401 .12040407 .26950403 .44300-00 
.33000400 .37570401 .11510407 .33160403 . 44650-00 
.33500400 .37010401 .11210407 .30340403 .44700-00 
.34000400 .36470401 .11030407 .43130403 .44650-00 
.34500400 .35940401 .10090407 .44920403 .44120-00 
.35000400 .35430401 .10660407 .48050403 .43300-00 
.35500400 .34930401 .10500407 .49000403 .42500-00 
.36000400 .34440401 .10030407 .51370403 .41600-00 
.36500400 .33970401 .80380406 .56130403 .41360-00 
.37000400 .33510401 .61940406 .60350403 .41370-00 
.37500400 .33070401 .44190406 .60940403 .42200-00 
.38000400 .32630401 .29790406 .60800403 .43890-00 
.38500400 .32210401 .20440406 .60980403 .45750-00 
.39000400 .31790401 .15320406 .62390403 .47420-00 
.39500400 .31390401 .12420406 .69120403 .49000-00 
.40000400 .31000401 .95960405 .84990403 .50400-00 
.40500400 .30620401 .74530405 .99280403 .51500-00 
.41000400 .30240401 .63890405 .10740404 .52560-00 
.41500400 .29880401 .59030405 .11050404 .53550-00 
.42000400 .29520401 .56190405 .11040404 .54590-00 
.42500400 .29180401 .53410405 .10870404 .55400-00 
.43000400 .28840401 . 4840D40S .10680404 .56140-00 
.43500400 .28510401 .41110405 .11000404 .56850-00 
.44000400 . 28180401 .33990405 .12160404 .57500-00 
.44500400 .27860401 .29010405 .13100404 .58060-00 
.45000400 .27550401 .26750405 .13880404 .58590-00 
.45500400 .27250401 .24540405 .14350404 .59120-00 
.46000400 .26960401 .22380405 .14520404 .59580-00 
. 46500400 . 26670401 . 20260405 . 14510404 . 59970-00 
.47000400 .26380401 .18190405 .14510404 .60350-00 
.47500400 .26100401 .16160405 .14700404 .60660-00 
. 48000400 . 25830401 . 14930405 . 15030404 . 60970-00 
.48500400 .25570401 .14190405 .14430404 .61250-00 
.49000400 .25310401 .13420405 .14350404 .61520-00 
.49500400 .25050401 .12650405 .14540404 .61760-00 
.50000400 .24800401 .12020405 .14510404 .62000-00 
. 5050D400 .24550401 .11430405 . 1440D404 .62200-00 
.51000400 .24310401 .10810405 .14170404 .62400-00 
. 51500400 . 24080401 . 1020D405 . 13850404 . 62620-00 
.52000400 . 2385D401 .96910404 . 1390D404 .62810-00 
.52500400 .23620401 .92370404 .14100404 .63000-00 
.53000400 .23400401 .88330404 .14070404 .63210-00 
.53500400 .23180401 .84830404 .13940404 .63380-00 
.54000400 .22960401 .81250404 .13720404 .63520-00 
.54500400 .22750401 .77040404 .13540404 .63720-00 
.55000400 .22540401 .72900404 .13370404 .63880-00 
.55500400 .22340401 .70100404 .13360404 .64000-00 
.56000400 .22140401 .67710404 .13190404 .64110-00 
.56500400 21950401 .65250404 .13300404 .64260-00 

.57000400 .21750401 .62550404 .13380404 .64370-00 
.57500400 .21560401 .59910404 .13470404 .64480-00 
.58000400 .21380401 .57790404 .13470404 . 6458D-00 
.58500400 .21200401 .55960404 .13470404 .64690-00 
.59000400 .21020401 .54160404 .13410404 .64770-00 
.59500400 .20840401 .52080404 .13290404 .64860-00 
.60000400 . 2067D401 .49990404 .13200404 .64970-00 
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Table 2 AMI solar spectral irradiance , P, silicon absorption 

coefficient* a, and silicon surface reflection coefficient 
R(A)» as a function of photon wave length, A. 88.92 mW/cm 
(part 2, 6000-12640A) 


Hu) 

. 60000+00 
. 60500+00 
. 61000+00 
. 62000+00 
. 63000+00 
. 64000+00 
.65000+00 
. 66000+00 
. 67000+00 
.68000+00 
. 60000+00 
. 69800+00 
. 70000+00 
. 71000+00 
. 72000+00 
. 72800+00 
. 73000+00 
. 74000+00 
. 75000+00 
. 76200+00 
. 77000+00 
. 78000+00 
. 79000+00 
. 80000+00 
. 80600+00 
. 82500+00 
. 83000+00 
. 83500+00 
. 84600+00 
. 86000+00 
. 87000+00 
. 87500+00 
. 88700+00 
. 90000+00 
. 90700+00 
. 91500+06 
. 92500+00 
. 93000+00 
. 94000+00 
. 95000+00 
. 95500+00 
. 96500+00 
. 97500+00 
. 98500+00 
. 10180+01 
. 10820+01 
. 10940+01 
. 10980+01 
. 11010+01 
. 11280+01 
. 11310+01 
. 11370+01 
. 11440+01 
. 11470+01 
. 11780+01 
. 11890+01 
. 11930+01 
. 12220+01 
. 12360+01 
. 12640+01 


E(eV) 
20670+01 
. 20500+01 
. 20330+01 
.20000+01 
. 19680+01 
. 19370+01 
. 19080+01 
. 18790+01 
.18510+01 
. 18230+01 
. 17970+01 
. 17760+01 
. 17710+01 
.17460+01 
. 17220+01 
. 17030+01 
. 16990+01 
. 16760+01 
. 16530+01 
. 16270+01 
. 16100+01 
. 15900+01 
. 15700+01 
. 15500+01 
. 15380+01 
. 1S03D+01 
. 14940+01 
. 14850+01 
. 14660+01 
. 14420+01 
. 14250+01 
. 14170+01 
. 13980+01 
. 13780+01 
. 13670+01 
. 13550+01 
. 13410+01 
. 13330+01 
. 13190+01 
. 13050+01 
. 12980+01 
. 12850+01 
. 12720+01 
. 12590+01 
. 12180+01 
. 11460+01 
. 11330+01 
. 1129D+01 
. 11260+01 
. 10990+01 
. 10960+01 
. 1090D+01 
. 10840+01 
. 10810+01 
. 10530+01 
. 1043D+01 
. 10390+01 
. 10150+01 
. 10030+01 
.98100+01 


a(CB~ 1 ) 

. 49990+04 
.47950+04 
. 46270+04 
. 43000+04 
. 39820+04 
. 36750+04 
. 33770+04 
.31300+04 
.29060+04 
. 27140+04 
. 25290+04 
. 23850+04 
. 23500+04 
. 21790+04 
. 20330+04 
. 19200+04 
. 18920+04 
. 17540+04 
. 16200+04 
.14640+04 
. 13620+04 
. 12590*04 
. 11780+04 
. 11000+04 
. 10440+04 
. 87440+03 
. 83580+03 
. 80000+03 
. 72280+03 
. 63060+03 
. 57090+03 
. 54160+03 
. 47380+03 
. 41320+03 
. 38130+03 
. 34550+03 
. 30530+03 
. 28660+03 
. 24970+03 
.21360+03 
. 19680+03 
. 16990+03 
. 14350+03 
. 11770+03 
. 52000+02 
. 65680+01 
. 47500+01 
. 42100+01 
. 38500+01 
. 16000+01 
. 14600+01 
.11500+01 
. 88000+00 
. 76000+00 
. 62000-01 
. 25000-01 
. 19000-01 
. 40000-02 
. 20000-02 
. 00000-00 


P(W/m 2 M) 

. 13200+04 
.13110+04 
. 13080+04 
. 1294P+04 
. 12810+04 
. 12720+04 
. 12570+04 
. 12440+04 
. 12270+04 
. 12100+04 
.11960+04 
. 10100+04 
.11750+04 
. 11570+04 
. 11350+04 
. 10030+04 
. 11180+04 
. 10950+04 
. 1077D+04 
. 79400+03 
. 10390+04 
. 10190+04 
. 10000+04 
. 98120+03 
. 87440+03 
. 93160+03 
. 9218D+03 
. 91240+03 
. 47620+03 
. 50640+03 
. 45380+03 
. 44920+03 
. 44860+03 
. 44890+03 
. 45520+03 
.46150+03 
. 27900+03 
. 22180+03 
. 31340+03 
. 29650+03 
.32110+03 
. 34440+03 
. 5769D+03 
. 54460+03 
.61750+03 
. 51290+03 
4641D+03 
. 50370+03 
. 50480+03 
. 13510+03 
. 15220+03 
.14310+03 
. 19120+03 
. 17450+03 
. 39930+03 
. 40220+03 
. 42400+03 
.39180+03 
. 39080+03 
. 32920+03 


l-R(X) 

. 64970-00 
.65030-00 
. 65120-00 
. 65290-00 
. 65430-00 
. 65590-00 
. 65710-00 
. 65870-00 
.66000-00 
.66070-00 
. 66200-00 
.66320-00 
. 66380-00 
. 66480-00 
. 66600-00 
. 66700-00 
.66720-00 
. 66790-00 
. 66900-00 
.67010-00 
. 67050-00 
. 67100-00 
. 67260-00 
. 67350-00 
. 67380-00 
. 67460-00 
. 67500-00 
. 67520-00 
. 67580-00 
. 67660-00 
. 67710-00 
. 6774D-00 
. 67800-00 
. 67840-00 
. 67880-00 
.67910-00 
. 67950-00 
. 67970-00 
. 68000-00 
. 68050-00 
. 68070-00 
.68100-00 
.68130-00 
. 68130-00 
. 68270-00 
. 68420-00 
. 68450-00 
. 68470-00 
. 68480-00 
. 68500-00 
. 68520-00 
. 68530-00 
. 68540-00 
. 68550-00 
. 68600-00 
.68610-00 
. 68620-00 
. 68660-00 
. 68690-00 
. 68730-00 


For reference purposes* the available photocurrent as a function of 
the total silicon solar cell thickness under the AMI condition Is computed. 

One pass is assumed, that Is* there Is complete transmission at the Hack 
Interface. The result is tabulated below and serves as a guide for a quick 
estimate of the maximum possible short circuit current for a given cell thickness. 
Higher values of short circuit current for thin cells than the tabulated 
values are possible if appreciable reflection occurs at the back interface* 
giving rise to multiple Internal reflections and absorptions. 

TABLE 3 

MAXIMUM AVAILABLE PHOTOCURRENT AS A FUNCTION OF THE THICKNESS 
OF THE SILICON SOLAR CELLS 


THICKNESS 

(m) 


AVAILABLE PHOTOCURRENT 
(raA/cm^) 


0.1 

0.2 

0.5 

1.0 

2.0 

5.0 

10.0 

20.0 

50.0 

100.0 

150.0 

200.0 

250.0 

300.0 

400.0 

500.0 
1000.0 
2000.0 
5000 

10000 


3.1318 

4.7755 

8.1773 

11.812 

16.219 

22.296 

26.208 

29.088 

31.491 
32.772 
33.438 
33.880 
34.198 
34.353 
34.826 
35.101 
35.943 
36.839 
37.702 

37.491 
40.853 
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The material properties required in the silicon solar cell performance 
calculations are the intrinsic carrier density and the mobilities of electrons 
and holes. A more sophisticated model would also require the energy levels 
of the shallow level dopant impurities (phosphorus and boron) in silicon as 
well as their dependences on their concentration. These Impurity deionisation 
and impurity banding effects are not Included since reliable theoretical models 
and experimental data are not available. Their inclusion will only slightly 
change the results. • 

The intrinsic carrier concentration is computed from 
n ± - 2.51xl0 19 (m*/m) 3/2 (T/300) 3/2 exp(-E G /2kT) 
where the effective mass and the energy gap are given by 
m*/m * (m^iOp/m 2 )^ 2 * 0.81577 

+ 3.4353xl0’ 3 [l - (T/437.65) + (T/814.2) 2 + (T/1356) 3 ]T 
E g - 1.2080 - 2.800xl0~ 4 + DE G (eV) 

DE g - - q 2 /4Tre(L D +a) « - 1. 2206374/ (Lp+a) (eV/A) 

Lp - i/ekT/q^ (N+P) 

a - [3/4ir(!H-P)] 1/3 
e - 11.8x8.85419xlO“ 14 F/cm 
q - 1.602192xl0 -19 C 
k - 8. 617087x10” 5 V/K 

Debye screening which reduces the energy gap is given by DEg above. This 
expression is consistent with that previously used by Morin and Malta and 
attributed to Herring [22] but differs from those by recent authors based on 
less physical derivations. Debye screening is the only mechanism which can 
give a true energy gap narrowing. Other mechanisms, such as the random position 
distribution of impurity and host atoms, can give apparant energy gap narrowing 
known as band tail states. However, these tail states are localized or partially 
localized low mobility or zero mobility states and they cannot be Included in 
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which ar« carriers In tha band with high mobilities. Energy gap narrowing 

la not included in this calculation, however, its effect is only Important in 

tha emitter region which has a concentration of majority carriers. Since 

emitter recombination is not the dominant factor (base recombination is) on 

the efficiency of sine doped sllcon solar cells, the inclusion of energy gap 

narrowing will reduce the performance only very slightly* This would not be 

the case if the base dopant impurity concentration is Increased substantially, 
17 -3 

for example, to 10 cm or above. 

The mobilities for majority and minority carriers are assumed to be 
equal and simple analytical fits are made to the experimental data. For AMI 
operation, the current density is sufficiently low that even a factor of two 
change of mobility would not significantly affect the computed results. The 
mobility fits given below, however, are better than 10%. The scattering by 
the lattice vibrations (acoustical and optical or intervalley phonons) and by 
ionised impurities are taken into account but the scattering by dislocations 
and neutral impurities are not included whose effects are expected to be small. 
The three scattering mechanisms are combined using the simple Matthiessen's 
rule 

"' 1 * “a 1 + “r 1 + ^ 

which give excellent fits to the published experimental data. More elaborate 
formulae based on complicated theoretical models do not justify the additional 
computation time required compared with this simple Mathiessen model. The 
three contributions are: 


2. 18x10 V 1 * 5 (cm 2 /V-s) 

Electrons 

1. 30x10 V 1 * 5 

Holes 

l^xlO 11 ^ 3 * 13 

Electrons 

6.64xl0 10 T“ 3 * 25 

Holes 

90(1 + (2.0xl0 18 /C I )(T/300)] 

Electrons 

45 [1 + (1.2x1018 /Ci ) (t / 3 oo)] 

Holes 
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In the Impurity scattering mobility formulae, the total ionized Impurity 
concentration, Cj, includes both the ionized shallow dopant impurities and 
the ionized deep level recombination centers, such as zinc in this case: 

Cj(x) - N^x) + ^(x) + HjC*) 

Where for zinc, N T (x) - N, + 2 M , where is the concentration of the singly 
charged zinc center and Nj is the concentration of the doubly charged zinc 
center. 

The numerical constants given in these mobility formulae were obtained 
by fitting the Mathiessen's rule to published experimental data of nearly pure 
silicon to heavily doped silicon (10^ to 10^® Impurity ions/cm^) in the 
temperature range of 4.2K to 600K. These fits are as good as those using 
more complicated theoretical formulae in recent publications. 
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